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We suggest a new type of the magnetic flux detector which can be optimized with respect to the 
measurement back-action, e.g. for the situation of quantum measurements. The detector is based 
on manipulation of ballistic motion of individual fluxons in a Josephson transmission line (JTL), 
with the output information contained in either probabilities of fluxon transmission/reflection, or 
time delay associated with the fluxon propagation through the JTL. We calculate the detector 
characteristics of the JTL and derive equations for conditional evolution of the measured system 
both in the transmission/reflection and the time-delay regimes. Combination of the quantum- 
limited detection with control over individual fluxons should make the JTL detector suitable for 
implementation of non-trivial quantum measurement strategies, including conditional measurements 
and feedback control schemes. 

PACS numbers: 



I. INTRODUCTION 



One direction of current efforts aimed at development 
of mesoscopic solid-state qubits is realization of effective 
schemes of qubit measurements. Different versions of the 
qubits which use quantum dynamics of magnetic flux in 
superconducting loops Q, S H S H B EM are at 
the moment among the most advanced solid-state qubits. 
Practically all of them employ the measurement scheme 
based on the modulation by the measured qubit of the 
decay rate of the supercurrent of a Josephson junction 
or a SQUID. This process can be viewed as tunneling of 
a magnetic flux quantum and has several attractive fea- 
tures as the basis for measurement. Most important one 
is sufficiently large sensitivity which comes from strong 
dependence of the tunneling amplitude on parameters of 
the tunneling potential controlled by the measured sys- 
tem. There is, however, an important disadvantage of 
this approach. In simple few-junction structures, the 
supercurrent decay brings the detector into the finite- 
voltage state characterized by large energy dissipation, 
which strongly perturbs both the system and the detec- 
tor itself, and makes it impossible to repeat the mea- 
surements sufficiently quickly. This strong perturbation 
also destroys the state of the measured quantum system 
by changing it in an uncontrolled way. Both of these 
factors prevent realization of non-trivial quantum mea- 
surement strategies which are typically based on contin- 
uous measurements or a sequence of successive discrete 
measurements. The goal of this work is to suggest and 
analyze a new flux detector which should not have this 
drawback. While still based on tunneling of individual 
magnetic flux quanta, the detector nevertheless avoids 
the transition into the dissipative state by using ballistic 
motion of the flux quanta |TjJ . The detector should have 
the quantum-limited back-action on the measured system 
and time resolution sufficient for realization of non-trivial 
quantum measurements of superconducting qubits. The 
cost of achieving this is the need for single-flux-quantum 
(SFQ) support electronics required for operation of 



this detector. Adaptation of the SFQ circuits to qubit 
applications 0, flH ITr| is an important, and not fully 
solved, problem of development of scalable superconduct- 
ing qubits. 



II. JOSEPHSON TRANSMISSION LINE AS 
FLUX DETECTOR 

The suggested detector is based on the ballistic mo- 
tion of fluxons in the Josephson transmission line (JTL) 
formed by unshunted junctions with critical currents Ic 
and capacitances C coupled by inductances L as shown in 
Fig. 1. The detector can be viewed as the flux analog of 
the quantum point contact (QPC) charge detectors (see 
[l7|| and references therein) used for measurements of the 
quantum-dot qubits. Both detectors make use of the abil- 
ity of the measured system to control the ballistic motion 
of independent particles (electrons in the QPC, fluxons 
in the JTL detector) through a one-dimensional channel. 
The JTL detector should, however, provide more control 
over the propagation of individual fluxons than is possible 
with electrons in the QPCs. This leads to new measure- 
ment regimes (e.g., time-delay measurement) impossible 
with the QPC. 

The JTL detector uses the fact that the flux & e \x) 
generated by the measured system creates potential U(x) 
for the fluxons moving in the JTL (Fig. 1). The fluxons 
are injected one by one, with the period 1//, into the 
JTL by the generator, and are scattered by the potential 
U{x) localized in some region of the JTL. The injection 
frequency / is sufficiently low so that only one fluxon at 
a time moves inside the JTL. The fluxon scattering char- 
acteristics: transmission probability or the time delay 
associated with the motion through the JTL, are regis- 
tered by the receiver. Since these characteristics depend 
on the potential U{x) controlled by the measured system, 
they contain information about the state of this system. 

Quantitatively, Hamiltonian of the uniform JTL can 
be expressed in terms of the electric charges Q n and the 
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FIG. 1: (a) Equivalent circuit of the flux detector based on the 
Josephson transmission line (JTL) and (b) diagram of scat- 
tering of the fluxon injected into the JTL with momentum 
k by the potential U(x) that is controlled by the measured 
qubit. The fluxons are periodically injected into the JTL by 
the generator and their scattering characteristics (transmis- 
sion and reflection coefficients t(k), r(k)) are registered by the 



phases <f> n of the Cooper-pair condensate at each of its 
nodes n: 

H = J2 [Q 2 J 2C + - cos ^™) 



(1) 



In this Hamiltonian, Ej = hlc/2e is the Josephson cou- 
pling energy of the junctions, where Ic is the junction 
critical current, El = ($o/27r) 2 /2L is the characteristic 
magnetic energy of inductances L, where $o = Tr^/e is 
the magnetic flux quantum, and <f>n = 27r<I>^ e ' ) /$o is the 
phase difference across the inductance L of the nth seg- 
ment induced by the external magnetic flux through 
this segment. The charges and phases at each node are 
the conjugate variables that satisfy canonical commuta- 
tion relations: [</>„, Q m ] = 2eid nm . 

In what follows, we are interested in the regime of small 
inductances, L <C <&q/Ic, in which the JTL acts as a uni- 
form ballistic channel for the fluxons. The phase differ- 
ence across each segment of the JTL is then small, and 
one can replace the phases (j) n with continuous function 
4>(x) of dimensionless coordinate x — n along the JTL. 
(Condition of validity of this approximation is discussed 
more quantitatively at the end of this Section.) In this 
regime, the JTL is equivalent to a distributed Josephson 
junction, and the Hamiltonian reduces to the stan- 
dard "sine-Gordon" Hamiltonian, which can be written 
as follows: 

H = 2E L J dx\^[{d T <pf + (d x <f>) 2 } + 

+Xj 2 (l-cos<f>)-^(x)d x 4>}. (2) 

Here t = ct is the time t normalized to the veloc- 
ity c = \j\jLC of propagation of excitations along 
the JTL in the absence of Josephson tunneling, and 
A,/ = (h/2eIcL) 1 / 2 is the Josephson penetration length 



of the junction. We note that in the notations used in 
Eq. (J2J) and below, all distances along the JTL, includ- 
ing x and Aj, are dimensionless, and are measured in 
units of the cell size a of the JTL (see Fig. 3). Conver- 
sion of most equations to absolute units of length can be 
achieved simply by changing interpretation of all quanti- 
ties (e.g., L, C, Ic) from "per cell" to "per unit length". 

The commutation relations for the charges and phases 
give the following equal-time commutation relations for 
the field 0(x,t) in the Hamiltonian ©: 



[</>{x),d T <l>{x')] = P 2 S(x-x' 



(3) 



where the parameter f3 2 = (4e 2 /h) \/L/C measures the 
wave resistance y/L/C of the JTL in the absence of 
Josephson tunneling relative to the quantum resistance. 
Known results for the quantum sine-Gordon model (see, 
e.g., 13) show that when f3 2 > 8tt, i.e., y/L/C > 
h/e 2 ~ 25 kfi, quantum fluctuations of the field 4> com- 
pletely destroy the quasiclassical excitations of the junc- 
tion. The transition at 1 — 8ir should be qualitatively 
similar to the analogous resistance-driven transition in 
small Josephson junctions 0. This analogy suggests 
that the dynamics of the supercurrent flow in the JTL 
with large wave resistance should be described in terms 
of tunneling of individual Cooper pairs |2(| . While this 
limit might be reachable in very narrow and thin JTLs 
of sub- micron width |2l| . we assume here more typical 
situation, when yf 'L/C is on the order of 10 - 100 Ohm 
and j3 2 <C 1. In this case, the JTL supports a number 
of quasiclassical excitations including, most importantly 
for this work, topological solitons that carry precisely one 
quantum of magnetic flux each. Dynamics of such "flux- 
ons" is e quiv alent to that of stable, in general relativistic, 
particles 18] with the terminal velocity c = 
mass 



m ~ &huj p /c 2 p 2 = (2h/e) 3 / 2 (I c L) 1 / 2 C , 



(4) 



where uj p — (2eIc/fiC) 1 ^ 2 — c/Xj is the plasma fre- 
quency. Another type of quasiclassical excitations in the 
JTL are the small-amplitude plasmon waves with fre- 
quency 



u{k) 



-c 2 k 2 )^ 2 . 



(5) 



for a wavevector k. 

In this work, we are interested in the "non-relativistic" 
regime of fluxon dynamics, when velocity u of its motion 
is small, u <§; c. Equations Q and JSJ show that in this 
regime, the fluxon kinetic energy e = h 2 k 2 /2m, where 
k = mu, can be made smaller than the lowest plasmon 
energy hu p , 



huj p ■ (2u/cP) 2 . 



(G) 



so that for u < c(3/2 the fluxon can not emit a plas- 
mon even when it is scattered by non-uniformities of 
the JTL potential J^]. Intrinsic dissipation associated 
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with emission of plasmons is then suppressed, and the 
fiuxon motion in the JTL should be elastic, provided that 
other, "extrinsic" , sources of dissipation are also suffi- 
ciently weak. Although the JTL operation as the flux 
detector should be possible for moderately strong fiuxon 
dissipation, the dissipation would prevent the detector 
from reaching the quantum-limited regime. 

The shape of the scattering potential U(x) for fluxons 
created by the measured system is determined by the 
convolution of the distribution of the flux ^>^(x) with 
the distribution of current in each fiuxon |23l l24j , and 
can be written as: 



U(x) 



2vrL 



where <j>o(x) is the shape of the fiuxon that in general 
can be distorted by the potential U(x) itself. If, how- 
ever, the potential is smaller than uj p , or does not vary 
appreciably on the scale of the size of the fiuxon given 
by Aj, the changes in the fiuxon shape are negligible, 
and one can use in Eq. (0 the regular fiuxon shape in 
the uniform case, which in the non-relativistic limit is 
4>o(x) = 4tan _1 [exp(a;/A,/)]. One of the implications of 
Eq. Q is that the width of the scattering potential U (x) 
can not be made smaller than Xj. 

If the measured system is coupled to the JTL not in- 
ductively (as in Fig. 1) but galvanically, and injects the 
current j^ (x) in the nodes of the JTL, potential created 
by this current for the fluxons is still given by Eq. J7J if 
one makes the substitution 



{l/L)d<S> [e \x)/dx= j (e) (x). 



(8) 



Discrete version of Eq. JSJ is illustrated in Fig. 2. Ex- 
ternal flux <E>! e ' > through the nth cell of the JTL induces 
circulating current j n = /i in this cell. The dif- 
ference between the currents j n in the neighboring cells 
creates the currents j n through the Josephson junctions 
of the JTL: jh ^ — j n —jn—x- I n the continuous limit, this 
relation gives Eq. (J8J. It implies that direct injection of 

(e) 

the currents j n in the JTL junctions indeed produces 
the same distribution of phases and currents as mag- 
netic coupling generating the fluxes <&n^- Throughout 
this work we assume that the measured system does not 
inject the net current in the JTL, J dxj^ e \x) = 0, i.e., 
the corresponding flux JSJ and the potential J7J vanish 
outside the coupling region: <J>( £ ) (x) -> for x -> ± 

Although the assumed condition 1 <g; 1 makes quan- 
tum fluctuations of the fiuxon shape small, the dynamics 
of the fiuxon as a whole can still be completely quan- 
tum. An example of quantum dynamics of fluxons of 
this type was recently observed experimentally 25] . JTL 
detector uses the quantum fiuxon dynamics in the po- 
tential described by Eqs. J7J) or JSJ. As was mentioned 
above, operation of the detector requires that the flux- 
ons are injected by the generator into one end of the JTL 
and observed by the receiver at the other end (Fig. 1). 




FIG. 2: Equivalence between the magnetic and galvanic cou- 
pling to the detector. External fluxes through the JTL 
cells induce the circulating currents j n in them which in their 
turn create the currents jn = j„ — j n -i through the JTL 
junctions. The situation would be the same if the currents 
jn 6 ' are injected directly into the JTL junctions by external 
system. 



Both of these circuits can be designed follo wing the gen- 
eral principles of the SFQ electronics 0, Il5j | and will 
not be discussed explicitly here. We simply assume that 
the ends of the JTL are matched appropriately to these 
circuits so that the fluxons can enter and leave the JTL 
without reflection and are injected in the JTL in an ap- 
propriate quantum state. This initial state ip(x, t = 0) of 
an injected fiuxon is characterized by the average fiuxon 
velocity u and the wave-packet ipo(x) defining its posi- 
tion: 



tfj(x, t = 0) = ipo(x)e lkoX , fc = mu. 



(9) 



As we will see from the discussion in the next Section, 
many properties of the JTL detector are independent of 
the specific shape of the wave-packet ipo(x), as long as 
it is well localized both in coordinate space and space of 
the fiuxon momentum k. They depend only on the wave- 
packet width £ in coordinate space and the corresponding 
width 5k ~ l/£ in momentum space. These parameters 
should satisfy two obvious conditions: £ <C I, where / is 
the total length of the JTL, and 5k <C feci. We assume a 
stronger form of the second condition that follows from 
the requirement that the broadening of the wave-packet 
by Sx ~ Skt/m because of the momentum uncertainty 5k 
during the typical time t ~ l/u of the fiuxon propagation 
through the JTL is negligible in comparison to the initial 
width: Sx <C £• The discussion below shows that this 
requirement is necessary for the quantum-limited opera- 
tion of the JTL detector in the time-delay mode. The two 
conditions mean that the width £ is limited as follows: 



(Z/fe ) 1/2 «£«Z. 



(10) 



In cases when it will be necessary to specify the shape 
of the wave-packet ^po{x) we will take it to be Gaussian: 



-i/4 e ~(x-x) 2 m 2 



(11) 



where x is the initial fiuxon position in the JTL. Be- 
sides being well localized as necessary both in the mo- 
mentum and coordinate space, the wave-packet l|l(J|) can 
be obtained as a result of the fiuxon generation process 
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that can be implemented naturally with the SFQ circuits. 
This process consists of the two steps: relaxation of the 
fluxon to the ground state of a weakly damped and nearly 
quadratic Josephson potential with the required width £ 
of the wavefunction of the ground state, and then rapid 
acceleration to velocity u after this potential is switched 
off. 

While the assumed condition /3 2 < 1 does not preclude 
the quantum dynamics of fluxons in the JTL, Eq. Q 
shows that decreasing increases the fluxon mass. This 
makes it more difficult to maintain quantum coherent 
dynamics of fluxons, which in the case of large mass be- 
comes more susceptible to perturbations. For instance, to 
avoid the effects of discreteness on the fluxon dynamics in 
the case of the discrete JTL structure (Fig. 1), the wave- 
length of its wavefunction © should be larger than the 
size of one cell. This condition can be written as fco < 1 
and means that the fluxon velocity and kinetic energy can 
not be larger than, respectively u = H/m = c ■ (/3 2 A,//8) 
and e = ti 2 /2m, i.e. 

e = hc[3 2 \j/16= (e 2 /C)(A, 7 /4). (12) 

We see that both the fluxon velocity u and its kinetic en- 
ergy e decrease with decreasing /3 2 and any external per- 
turbation (e.g., fluctuations of the critical current of the 
JTL junctions) affects the fluxon motion more strongly 
at small 0. For realistic values of Aj and (3 the limita- 
tion on the fluxon energy l|12|) is stronger than © . Since 
the fluxon size Aj can not be increased much for obvi- 
ous reasons of convenience of fluxon manipulation, the 
energy |Q depends only on the junction capacitance C, 
and can be increased only by decreasing the junction size 
(see the discussion below). In practical terms, the fluxon 
energy (|12[1 should be at least larger than temperature 
for the fluxons to be generated by an SFQ generator. All 
this means that the JTL can operate as the quantum de- 
tector only with junctions of sub-micron size. Although 
the estimate (|12f) was made for the discrete JTL, the 
same practical conclusion would be reached in the case 
of nominally-uniform long Josephson junctions, for which 
the limitation on the fluxon momentum fco would be set 
by unavoidable spatial fluctuations of the junction pa- 
rameters. 

We end this Section with a discussion of experimental 
realization of the JTL detector with conventional thin- 
film "tri-layer" fabrication technology of the SFQ cir- 
cuits. Simplified but sufficiently realistic geometry of one 
JTL cell in this case is shown in Fig. 3. It contains 4 su- 
perconductor (niobium) films, the lower and the top of 
which (ground plane and the counter-electrode) are used 
to define inductances, while the two layers in the middle 
are separated by a very thin layer of insulator (the tri- 
layer structure) and are used to produce the Josephson 
junctions. Such a discrete multi-layer configuration of the 
JTL is more convenient than a uniform long Josephson 
junction, since it enables one to optimize separately the 
junctions and inductances of the JTL, most importantly, 
to reduce effective capacitance. 
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FIG. 3: Geometry of one cell of the multi-layer JTL: top view 
and vertical cross-section. Internal layers of the structure 
are used to create Josephson junctions (JJs), while the two 
external layers, the ground plane and the counter-electrode, 
separated vertically by the distance ti, define inductance L. 
The length a of the cell is set by the distance between the 
junctions. 

One possible drawback of the discrete configuration 
could be the parasitic periodic potential u(z) for the flux- 
ons produced by the JTL discreteness (see, e.g., |26j). 
This potential, however, is negligible when Aj is not too 
small, Aj > 2. Indeed, using the Hamilonian Q and 
equation for the static distribution <p n that follows from 
it, we can write the energy E of the static fluxon as 

E = E J '^2e n , e n = 1 - cos<p n - -0„sin0 n . (13) 

n 

Applying the Poisson summation formula to this expres- 
sion, we get in the limit of sufficiently large A,/: 

E = E [l + acos(27rz)] , (14) 

where Eq = mc 2 is the fluxon rest energy, z is the fluxon 
coordinate in units of the JTL period a, and 

a= l -J dye( y y 2 *y x > . (15) 

Here e(y) is a smooth envelope of the discrete distribution 
e n (|13f) associated with the static fluxon <f> n . Since e(y) is 
smooth in the limit of large Aj, Eq. (|15fl implies that the 
amplitude Se of the z-dependent part of the fluxon energy 
is exponentially small in Aj. Specifically, if e(y) is taken 
as for the regular "continuous" fluxon <j)(y) = 4tan _1 [e y ], 
we see that 

a cx e^ Xj , 

and periodic potential associated with the JTL discrete- 
ness is completely negligible for A,/ > 2 even on the scale 
of small fluxon kinetic energy l|12|) . For operation of the 
JTL detector it is obviously convenient to keep Aj (i.e., 
the fluxon size) as close as possible to a minimum con- 
sistent with this inequality. 

In the multi-layer fabrication technology (Fig. 3), the 
Josephson junctions have the critical current density j c 
that can be varied within wide range, 10 -r 10 4 A/cm 2 , 
while specific junction capacitance cj changes only a lit- 
tle, from 30 fF //im 2 to 60 fF //an 2 , in this range of j c 's. 
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This means that there is a strong limitations on the min- 
imal junction capacitance C set by the junction diameter 
D, e.g., C ~ 30 fF for D ~ 1.0 /im. Inductances in the 
multi-layer JTL (Fig. 3) are characterized by the specific 
inductance per square K = fMy ■ (t j + 2A), where A is the 
superconductor penetration depth, and parasitic capaci- 
tance cl ~ eeoAi per unit area. For realistic « 0.2-^0.4 
/im, k and cl are 0.5 -j- 0.7 pH and 0.2 0.1 fF//um 2 , re- 
spectively. 

The discussion leading to Eq. (|12|) shows that the main 
requirement on the JTL parameters optimizing its opera- 
tion as the quantum-limited detector consists in making 
the junction capacitance C as small as possible. After 
this, the inductance L can also be reduced as long as this 
reduction is at least partly compensated by increase in j c 
and does not make the characteristic fluxon size Aj too 
large. The actual fluxon size is about 4 Aj, and since it is 
impractical to have significantly more than about 30 cells 
in the JTL, Aj is basically fixed within the range from 2 
to 3. For relatively small junction size D = 0.3 [im and 
high critical current density j c = 10 kA/cm 2 , the value 
of Aj within this range corresponds to L ~ 4 pH, induc- 
tance that is obtained for the width of the superconduc- 
tor strip w = 1 /im, and the cell size a ~ 5 fixa. For these 
parameters, the parasitic capacitance cl gives noticeable 
contribution (about 30%) to the effective junction capac- 
itance C ~ 4 fF. The propagation velocity along such a 
JTL is c ~ 10 13 cells/s, i.e., 5 ■ 10 7 m/s, and the wave re- 
sistance is about 30 Ohm giving 1 ~ 0.03. Useful fluxon 
velocity in the quantum-limited regime [estimated from 
Eq. is 10 11 cells/s, so that the total time of the 

fluxon propagation through the JTL (which limits the 
detector time resolution) is about 0.3 ns. We se that the 
JTL detector should have very good time resolution even 
in the quantum-limited regime, but requires in this case 
the JTL parameters which are close to the limit of the 
current SFQ fabrication technology. 

If one abandons the goal of realizing the quantum- 
limited detection, requirements on the JTL parameters 
become much more routine. A typical set of parameters 
D ~ 1.5 /im, j c = 1 kA/cm 2 , and w — 3/im, a = 10 /im, 
produces the JTL with c ~ 3T0 12 cells/s, wave resistance 
about 2 Ohm, and Aj ~ 2. This should give the detec- 
tor that is not quantum-limited but is very fast and has 
the time resolution on the order of 0.1 ns. As discussed 
in Sec. 5, such a detector could be used in non-trivial 
schemes of quantum measurement despite classical dy- 
namics of fluxons in it. 



III. MEASUREMENT DYNAMICS OF THE JTL 
DETECTOR 

As discussed above, the measurement by the JTL de- 
tector consists qualitatively in scattering of the fluxons 
in the JTL by the potential controlled by the measured 
system. This process has the simplest dynamics if the 
measured system is stationary. In this case it is conve- 



nient to consider it in the basis of the eigenstates \j) of 
the system operator (e.g., magnetic flux in the qubit loop 
in the example shown in Fig. 1) which couples it to the 
JTL. In each state \j), the system creates different poten- 
tial Uj(x) for the fluxons propagating through the JTL. 
Different realizations Uj (x) of the JTL potential produce 
different scattering coefficients for injected fluxons: the 
amplitude tj(k) of transmission through the JTL, and 
the amplitude rj(k) for the fluxon to be reflected back 
into the generator. Both amplitudes are in general func- 
tions of the fluxon momentum k. Since they depend on 
the state \j) of the measured system, scattered fluxons 
carry information about 

In general, the process of quantum measurement can 
be understood as creation of entangled state between the 
measured system and detector as a result of interaction 
between them. The states of the detector are quasiclassi- 
cal and suppress quantum superposition of different out- 
comes of measurement. The two consequencies of this 
process are the acquisition of information about the sys- 
tem by the detector and "back-action" dephasing of the 
measured system - see e.g., [2?} . For the JTL detec- 
tor, the detector-system entanglement arises as a result 
of fluxon scattering, and the rates of information acquisi- 
tion and back-action dephasing can be expressed in terms 
of the scattering parameters. In this respect, the JTL is 
very similar to the QPC detector and its characteristics 
can be obtained following the derivation [13 for the QPC. 

Evolution of the density matrix p of the measured sys- 
tem in scattering of one fluxon can be obtained by consid- 
ering first the time dependence of the total wavefunction 
of the fluxon injected into the JTL and the wavefunction 
J2j c j\j) °f the measured system: 

1>{x, t = 0) ■ c i\j) - E °i m*, *) • U) ■ ( 16 ) 

i j 

Here the initial fluxon wavefunction ip(x, t = 0) is given 
by Eq. @ and its time evolution ipj(x, t) depends on the 
realization Uj(x) of the scattering potential created by 
the measured system. Evolution of p is obtained then by 
tracing out the fluxon part of the wavefunctions (|16fl : 

Pij = c t c* -> CiC* J dxipi(x,t)ip*(x,t) . (17) 

Qualitatively, the time evolution in Eq. (|16|l de- 
scribes propagation of the initial wavepacket © to- 
wards the scattering potential and then separation of 
this wavepacket in coordinate space into the transmit- 
ted and reflected parts that are well-localized on the op- 
posite sides of the scattering region. If we assume that 
the scattering potential U(x) has simple shape (for in- 
stance, does not have narrow quasi-bound states) and is 
non- vanishing only in some small region of size rj <C I, 
the time t sc from the fluxon injection to completion of 
the scattering process is not drastically different from 
the time l/u of free fluxon propagation through the JTL. 
Then at time t > t sc , the separated wavepackets move in 



6 



the region free from the j-dependent scattering potential 
and the unitarity of the quantum-mechanical evolution 
of ipj(x,t) implies that the overlap of the fluxon wave- 
functions in Eq. (|17|l becomes independent of t. 

This overlap can be directly found in the momentum 
representation: 



dxipi(x, t)ip*(x, t) 



dk\b{k)\ 2 [m+r ir *]. (18) 



Here b(k) is the probability amplitude for the fluxon to 
have momentum k in the initial state 10, e.g., in the case 
of the Gaussian wavepacket iJTTJ, 



b(k) = (fA) 



-1/4 (k-k ) 2 e/2-i(k-k )x 



(19) 



Equations i|17fl and (|18fl show that the diagonal ele- 
ments of the density matrix p do not change in the pro- 
cess of scattering of one fluxon, while the off-diagonal 
elements are suppressed by the factor 



dklbik^Uk^i^ + niky^k)} 



< 1 



(20) 



The inequality in this relation can be proven as the 
Swartz inequality for the scalar product in the Hilbert 
space of "vectors" {tj(k), rj(k)} weighted by |6(fc)| 2 . The 
suppression of the off-diagonal elements of p due to the in- 
teraction with fluxons is manifestation of the back-action 
dephasing of the measured system by the JTL detector. 
If we add the suppression factors for the fluxons injected 
into the JTL with frequency /, the rate of this dephasing 
is obtained as 



-/In 



dk\b(k)\ 2 [ti(k)t*{k) + r t (k)r*(k)} 



(21) 



Equation (|21|l is similar, but not identical, to the back- 
action dephasing rate by the QPC detector |l7|. The 
main difference with the QPC is the stage at which the 
summation over the momentum k is carried out. This 
difference reflects the fact that in contrast to electrons in 
the QPC, different /c-components of the fluxon wavefunc- 
tion do not scatter independently. They are constrained 
by the condition that one fluxon as a whole is either 
transmitted or reflected by the potential, since scatter- 
ing events of different fluxons are well-separated in time. 
As we will see below, this difference makes it possible to 
operate the JTL detector in the new, time-delay mode 
that is not possible with the dc-biased QPC detector. 

As examples of application of Eq. (|21|l we consider sev- 
eral specific cases motivated by the measurement regimes 
discussed below. In one, the phases of the scattering 
amplitudes are assumed to cancel out from Eq. (|21|l . 
while the variation of the absolute values of the am- 
plitudes with index j is small. The dephasing rate 
(|21|l can be expressed then in terms of the variations 
STj(k) of the fluxon transmission probability in differ- 
ent states \j) around some average transmission T(k): 



\t 3 {k)\ 2 = 5Tj{k) + T(k), 6Tj(k) < T{k). In the lowest 
non- vanishing order in 5Tj(k) we get: 



r y = / / dfc|6(fe)| 



8T(k)(l-T(k)) 



(22) 



This equation describes the "linear-response" regime of 
operation of the JTL detector, when its properties follow 
from the general theory of linear detectors ■ In par- 
ticular, the dephasing (|22|l can be understood as being 
caused by the back-action noise arising from the random- 
ness of the fluxon transmission/relection process. 

In the "tunnel" limit of weak transmission, |fj(fc)| <C 
1, and again assuming that the phases of the scattering 
amplitudes cancel out, Eq. 121fl reduces to: 



r« = if/2) 



dk\b{k)\ 2 {Uk)\~Uk)\f 



(23) 



Both Eqs. \\T2l and 1)2311 have direct analogues in the case 
of the QPC detector7|2lE!II3. 

As the last example that does not have the analogue 
in the QPC physics, we consider the situation when the 
reflection amplitudes are negligible, rj(k) = 0, and the 
system- JTL interaction modifies only the phases Xj(fc) of 
the transmission amplitudes: tj(k) — e lXj ( k \ so that 



r« = -/in 



dk\b(k)\ 



2J[xi(k)-Xj(k)] 



(24) 



In coordinate representation, this means that the scat- 
tering potential affects only the position and, in general, 
the shape of the transmitted wavepacket. If the width 8k 
of the initial fluxon state in the momentum representa- 
tion is sufficiently narrow, and the phases Xj(k) do not 
vary strongly over this momentum range, they can be 
approximated as 

Xj( k ) = X]( k o) - (k - k^Xj , Xj = -x'j{k Q ). (25) 

One can see directly that this approximation neglects dis- 
tortion of the fluxon wavepacket in the scattering process, 
while taking into account its shift Xj along the coordinate 
axis. This shift can be directly related to the "time de- 
lay" Tj due to scattering: tj = Xj/u (which can be both 
negative and positive depending on the scattering poten- 
tial). The dephasing rate (|24() can then be conveniently 
written down in the coordinate representation in terms 
of the initial fluxon wavepacket ipo(x): 



Fij = -/in 



dxipo[x — Xi)ipQ(x — Xj 



(26) 



Equation Q26JI shows explicitly that the back-action de- 
phasing by the JTL detector arises from the entangle- 
ment between the measured system and the scattered 
fluxons which are shifted in time by an interval dependent 
on the state of the system. The degree of suppression of 
coherence between the different states of the measured 



7 



system is determined then by the magnitude of the rel- 
ative shift of the fluxon in these states on the scale of 
the wavepacket width. For instance, if the initial fluxon 
wavepacket is Gaussian ffTTf) . Eq. gives: 

r y =f{x i -x j f/Ae. (27) 

Back-action dephasing represents only part of the mea- 
surement process. The other part is information acqui- 
sition by the detector about the state of the measured 
system. In the case of the JTL detector, this information 
is contained in the scattering characteristics of fluxons, 
and the rate of its acquisition depends on specific char- 
acteristics recorded by the fluxon receiver. There are at 
least two different possibilities in this respect. One is 
to detect the probability of fluxon transmission through 
the scattering region (or, equivalently, the corresponding 
probability of the fluxon reflection back into the gener- 
ator, see Fig. 1). Another possible detection scheme is 
for the receiver to measure the time delay associated with 
the fluxon propagation through the JTL. Even if the mea- 
sured system changes the potential Uj (x) in such a way 
that the fluxon transmission probability is not affected, 
potential can still change the fluxon propagation time, 
which will contain then information about the state of 
the system. In general, one can have a situation when 
the information is contained both in the changes of the 
propagation time and transmission probability, and in or- 
der not to loose any information one would need to detect 
both scattering characteristics. In this work, we consider 
only the two "pure" cases of transmission and time-delay 
detection modes. 



A. Transmission detection mode 

If the detector records only the fact of the fluxon ar- 
rival at the receiver, then only the modulation of the 
fluxon transmission probability by the measured system 
conveys information about the system. The information 
contained in all other features of the scattering ampli- 
tudes (e.g., their phases or propagation time) is lost in 
the receiver. In this case, the rate of information ac- 
quisition can be calculated simply by starting with the 
probabilities of the fluxon transmission/reflection Tj and 
Rj, when the measured system is in the state \j): 

T 3 = J dk\b(k)\ 2 \t 3 (k)\ 2 , Rj = l-Tj. (28) 

Since the outcomes of successive fluxon scattering 
event are independent, the probability p{n) to have n 
out of N incident fluxons transmitted, is given by the 
binomial distribution Pj (n) = C^T^Rf- 71 , The task of 
distinguishing different states \ j) of the measured system 
is transformed into distinguishing the probability distri- 
butions Pj(n) for different js. Since the number N = ft 
of scattering attempts increases with time t, the distri- 
butions pj (n) become peaked successively more strongly 



around the corresponding average numbers TjN of trans- 
mitted fluxons. The states with different probabilities 
Tj can be distinguished then with increasing certainty. 
The rate of increase of this certainty can be character- 
ized quantitatively by some measure of the overlap of 
the distributions pj (n) . While in general there are differ- 
ent ways to characterize the overlap of different probabil- 
ity distributions |3l| . the characteristic which is appro- 
priate in the quantum measurement context [32L l33| is 
closely related to "fidelity" in quantum information |3l| : 
Zinb^fet")] 1 '' 2 - The rate of information acquisition 
(increase of confidence level in distinguishing states |i) 
and |j)) can then be defined naturally as |l7| : 

Wij = lnJ2\Pi(n)Pj(n)} 1/2 . (29) 

n 

Using the binomial distribution in this expression we get: 

= -fhxHTiTj) 1 ' 2 + [R t R 3 ) 1/2 ] , (30) 

where the transmission and reflection probabilities are 
given by Eq. J2HJ ■ 

Equation (|3(J[1 characterizes the information acqui- 
sition rate of the JTL detector in the transmission- 
detection mode. For an arbitrary detector, the infor- 
mation rate should be smaller than or equal to the back- 
action dephasing rate, and regime when the two rates 
are equal is "quantum-limited" . Comparing Eqs. I|21|) 
and H30fl one can see that for the JTL detector, indeed, 

< Tij , (31) 

and equality holds if several conditions are satisfied. First 
two conditions require that there is no information in the 
phases of the transmission amplitudes: 

fa{k) = M*>), Mk) = Bxg(tj(k)/ rj (k)), (32) 
Xj(k) - Xi(fc) = const. (33) 

The two conditions, ll-'il'l) and (|33|l . have different physical 
origin. Condition (|32|l implies that the scattered states 
contain no information on j that can be used in prin- 
ciple by arranging interference between the transmitted 
and reflected parts of the wavefunction [l^ ■ In practi- 
cal terms, the simplest way to satisfy this condition is to 
make the scattering potential symmetric Uj{—x) = Uj(x) 
for all states \ j). The unitarity of the scattering matrix 
for the fluxon scattering in the JTL implies in this case 
that <pj = tt/2 for any j. Condition l|33|l means that no 
information on j is contained in the shape and position of 
transmitted wavepackets that would be lost in the fluxon 
receiver operating in the transmission-detection mode. 
(Similar condition for the reflected wavepacket follows 
from Eqs. (|32() and H33[) .~) In general, condition l|33() re- 
quires that the spread 5k of the initial fluxon state over 
momentum gives rise to the uncertainty in the fluxon 
energy Se ~ u5k that is much smaller than the energy 
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scale ft of the transparency variation of the scattering 
potential Uj(x). 

One more condition of the quantum-limited operation 
is that the fluxon transmission probabilities are effec- 
tively momentum- and energy-independent in the rele- 
vant momentum range: 

\t 3 {k)\ 2 =Tj. (34) 

This condition again requires that 6e <C ft- It is more 
restrictive than the corresponding condition for the QPC 
detector which can be quantum-limited even in the case 
of the energy-dependent transmission probability 0,0, 
l35l |. To obtain Eq. (|34ll . one starts from the back-action 
dephasing rate l|21|) which can be written as 

r« = -/in J dk\b{k)\ 2 [\u{k)t 3 (k)\ + \n{k) rj (k)\] 

under the conditions ll-il'l) and . Swartz inequality for 
the functions |tj(fc)|: 

[TiTj] 1 / 2 > J dk\b(k)\ 2 \ti(k)tj(k)\ , 

and similar inequality for show that this Ljj and 

the information acquisition rate Wij satisfy the inequality 
(|31|l . Equality in Eq. could be reached only when 

\tj{k)\ =X j t(k), \r j (k)\ = X' j r(k), (35) 

where A's are some constants, so that the ratios 
and 1^(^)1/1^(^)1 are independent of k. 
Similarly to Eq. I|33(l , Eq. I|35|) demands that no informa- 
tion about the state of the measured system is contained 
in the shape of transmitted or reflected wavepackets. In 
general, when both the transmission and reflection prob- 
abilities are not small, the two relations in (|35|l are in- 
compatible. They have only a trivial solution in which 
all amplitudes are independent of k in the relevant range 
of fc's, thus proving Eq. ® for Tj ~ Rj ~ 1/2. The 
transmission and reflection probabilities have roughly the 
same magnitude when the fluxon energy e is close to the 
maximum U of the scattering potential. In this case, 
small spread of e: Se <C ft implies that the range r\ of the 
scattering potential U(x) should be small: r] <C £. 

Condition i|34|) of the quantum- limited operation of the 
JTL detector is not necessary when either Tj < 1 or 
Rj <C 1. In this case, one of the relations in l|35|) reduces 
to a trivial statement |fj (fc)| ~ 1 or |t/(fc)| — 1. The other 
relation gives then the actual condition of the quantum- 
limited operation that can be satisfied in principle with 
an arbitrary fc-dependent function t{k) or r(k). 

Finally, we show how Eq. (|34f> applies in the linear- 
response regime, when the variations STj of JTL trans- 
parencies between the different states \j) are small and 
back-action dephasing is given by Eq. (E2J). Expand mg 
Eq. (|3T)jl in 5T 3 : T 3 = T + 5T 3 , where all transparencies 
are defines as in Eq. we get: 

Wij = f(STi - STjf/[8T(l - Tj)] . (36) 



This equation differs from Eq. (|22|l only by the order in 
which the integration over momentum is performed. This 
means that the information and dephasing rates satisfy 
the inequality (|3 II) and are equal only if the transparency 
is constant in the relevant momentum range. 

In the linear-response regime, each individual fluxon 
carries only small amount of information, and it is conve- 
nient to employ the quasi-continuous description in which 
the fluxon receiver acts as the voltmeter registering not 
the individual fluxons, but the rate of arrival of many 
fluxons, i.e., the voltage V(t) across the junctions of the 
JTL. The average voltage in the state \j) is 

(V(t)) = /T,$o , (37) 

where (...) implies average over the scattering outcomes 
and over time t within the fluxon injection cycle. Be- 
cause of the randomness of the fluxon scattering, the ac- 
tual voltage fluctuates around the average values l|3"7|) 
even at low frequencies. The voltage fluctuations can be 
described as the shot noise of fluxons and its spectral 
density 

Sv(u) = J dre-^ T ((V(t + r)V(t)) - (V(t)) 2 ) (38) 

is constant at frequencies to below the fluxon injection 
frequency /. Straightforward calculation similar to that 
for the regular shot noise shows that this constant is 

So = fT(l - T)$g t (39) 

where in the linear-response regime we can neglect small 
differences STj of transparencies between the different 
states \j) in the expression for noise. This equation shows 
that in accordance with the general theory of linear quan- 
tum measurements, the information rate l|36|l can be in- 
terpreted as the rate with which one can distinguish dc 
voltage values (|37l) in thepresence of white noise with 
the spectral density l|3"§|) |27ll33l|. 

To summarize this subsection, we see that in the 
most relevant regime Tj ~ 1/2 which maximizes the 
detector response to the input signal, conditions of the 
quantum-limited operation of the JTL detector are given 
by Eqs. ll-il'l) and Q34JI. These conditions are satisfied if 
the scattering potential for the fluxons created by the in- 
put signal is symmetric and has the range r\ smaller that 
the size £ of the fluxon wavepacket. 

B. Time-delay detection mode 

Since the range r\ of the scattering potential Q) can 
not be smaller than the fluxon size Aj, it might be dif- 
ficult in practice to realize the condition r\ <^ £ needed 
for the quantum-limited operation of the JTL detector 
in the transmission-detection mode. For "quasiclassical" 
potential barriers that are smooth on the scale of fluxon 
wavepacket, 77 > £, the "transition" region of energies 
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near the top of the barrier where the reflection and trans- 
mission amplitudes have comparable magnitude, is nar- 
row. If the interval Se of the fluxon energies avoids this 
narrow region, then either transmission or reflection co- 
efficient can be neglected. Ballistic motion of fluxons in 
this regime still contains information about the potential 
Uj(x) that can be used for measurement. This infor- 
mation is contained in the time shift Tj caused by the 
propagation through the region of non-vanishing poten- 
tial. Quantum-mechanically, the time-shift information 
is contained in the phases of the scattering amplitudes. 
To be specific, we discuss here the regime when the JTL 
detector is operated in this "time-delay" detection mode 
using the transmitted fluxons, i.e., = 1. In the 

energy range where = 1, the same detection pro- 

cess is possible using the reflected fluxons, the advantage 
of the transmission case being the possibility to make 
use of the full range of values of the scattering potential: 
Uj(x) < and Uj(x) > 0. 

For sufficiently smooth potential Uj(x), the phase 
Xj(fc) of the transmission amplitude can be calculated 
in the quasiclassical approximation: 



can write the information acquisition rate of the JTL de- 
tector in the time-delay mode as follows: 



Xj (k)= / dx[2m{e-Uj{x))) x ' 2 



(40) 



If the potential is weak, Uj(x) <C e, potential-induced 
contribution to the phase 140|l is 



Xj (k) =-- I dxlfj(x) . 



(41) 



Under the adopted assumption of quasiclassical poten- 
tial and vanishing reflection, condition l|10|) of the negligi- 
ble broadening of the fluxon wavepacket still works in the 
presence of potential. In this case, one can use approxi- 
mation (|25|l for the phases Xji x ) which implies the shift 
of the wavepacket as a whole without distortion. The 
potential-induced part of the shift Xj — — x'j(ko) follows 
from Eq. I|40[) and has the classical form: 



Xj = / dx[l - — ^— ] , Uj {x) = [2(e - Uj{x))/m 
For weak potential, Uj(x) <C e: 



Xj ■. = -— I dxUj(x) 



1/2 _ 

(42) 
(43) 



Back-action dephasing rate by the JTL detector in this 
regime is given by Eqs. (|26(l and (|27|l . The information 
about the states \ j) contained in the shift Tj of the fluxon 
in time or, equivalently, coordinate Xj — TjU, can be read 
off by distinguishing different shifts Xj against the back- 
ground of the finite width £ of the fluxon wavepacket 
ipo(x)- Since l^oC^) | 2 gives the probability of finding 
the fluxon at coordinate x, this task is equivalent to 
the task of distinguishing two shifted probability dis- 
tributions (see Fig. 4) that was discussed above for the 
transmission-detection mode. Similarly to Eq. (|29|) . we 



Wo. 



/In / dx\if} (x - Xi)ip (x - Xj)\ . 



(44) 



Comparing this to the dephasing rate (|26|l . we see that in 
general the two rates satisfy the inequality l|31|l as they 
should. The rates are equal if the phase of the initial 
wavepacket tpo(x) of the injected fluxon is independent 
of x, i.e., if tpo(x) is essentially real. In particular, in the 
case of the Gaussian wavepacket (|19|l , the JTL detector is 
quantum-limited, Wij = , and the two rates are given 
by Eq. I|27l) . These considerations also imply that the 
JTL detector in the time-delay mode would lose the prop- 
erty of being quantum-limited if the fluxon wavepacket 
spreads noticeably during propagation through the JTL. 
This process creates non-trivial x-dependent phase of the 
wavepacket and makes the information acquisition rate 
smaller than the back-action dephasing rate. 



IV. CONDITIONAL EVOLUTION 

As is the case with any dephasing, the back-action 
dephasing by the JTL detector can be viewed as the loss 
of information. In the regime of the quantum-limited 
detection, the overall evolution of the detector and the 
measured system is quantum- coherent and the only 
source of the information loss is averaging over the 
detector. For a detector, different outcomes of the 
evolution are, however, classically distinguishable, and 
it is meaningful to ask how the measured system evolves 
for a given detector output. In the quantum- limited 
regime, specifying definite detector output eliminates 
all losses of information, and as a result there is no 
back-action dephasing present in the dynamics of the 
measured system conditioned on specific detector output. 




FIG. 4: Illustration of the information acquisition process 
by the JTL detector in the time-delay mode. The fluxon 
wave-packet is shifted by the distance Xj dependent on the 
state \j) of the measured system. In conditional description, 
observation of the fluxon with position x' changes the system 
wavefunction according to Eq. 



Conditional des crip tion in the quantitative form is ob- 
tained (see, e.g., [33, HE |23 ) by separating in the to- 
tal wavefunction the terms that correspond to a specific 
classical outcome of measurement and renormalizing this 
part of the wavefunction so that it corresponds to the 
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total probability of 1 . In the case of the JTL detector in 
the transmission-detection mode, for each injected fluxon 
there are two classically different outcomes of scatter- 
ing: transmission and reflection of the fluxon. Accord- 
ingly, the wavefunction of the measured system should 
be conditioned on the observation of either transmitted 
or reflected fluxon in each cycle of fluxon injection. The 
evolution of the total wavefunction "detector+measured 
system" during the scattering of one fluxon is described 
by Eq. (|To)) . If conditions and of the quantum- 
limited detection are satisfied, the transmission and re- 
flection amplitudes are effectively momentum indepen- 
dent [33 ■ Coordinate dependence of the scattered fluxon 
wavepakets is then the same in different states and 
can be factored out from the total wavefunction. The 
evolution of the measured system can then be condi- 
tioned on the transmission/reflection of a fluxon simply 
by keeping in Eq. I|16|) the terms that correspond to the 
actual outcome of scattering. If the fluxon is transmit- 
ted through the scattered region or reflected from it in a 
given injection cycle, amplitudes Cj for the system to be 
in the state \j) change, respectively, as follows: 

^T72> ^T72- ( 45 ) 

It is important to stress that the changes in the coef- 
ficients Cj for a system with vanishing Hamiltonian (as 
we assumed from the very beginning) is unusual from 
the point of view of Schrodinger equation, and provides 
quantitative expression of reduction of the wavefunction 
in the measurement process. Equations (|45|l are similar 
to those obtained in the "bayesian" approach 33] for the 
QPC detector in the tunnel limit. 

If dynamics of the JTL detector is not quantum- 
limited, then information is lost and dephasing is non- 
vanishing even in conditional evolution. To generalize 
Eqs. (|45|l to this case of finite "residual" dephasing, we 
need to start with Eq. i|17|) for the change of the density 
matrix pij of the measured system due to the scattering 
of one fluxon. To condition this change on the specific 
outcome of scattering, we limit the trace over the fluxon 
coordinate to the range containing only transmitted or 
reflected fluxons. In this way we see that /ay changes as 

Pij - Pij J rffc|6(fc)| 2 (t i (A ; )t*(fc))/E i p ii T j , (46) 

Pij -> Pij J &|6(fc)| 2 (r l (fc)r*(fc))/I] J p w i? J , (47) 

if the fluxon is, respectively, transmitted or reflected in 
a given cycle. Using the discussion of the transmission- 
detection mode in Sec. Ill, one can see immediately that 
the requirements H33fl and (|34|l of the JTL detector being 
quantum-limited, are equivalent to the density matrix py 
remaining pure in the conditional evolution described by 
Eqs. gSJ and g7J). 

The above derivation of conditional evolution equa- 
tions can be repeated with only minor modifications in 



the time-delay mode of the detector operation. In this 
case, different classical outcomes of measurements are the 
observed instances of time when the fluxon reaches the 
receiver, that for convenience can be directly translated 
into different fluxon positions x at some fixed time. If the 
fluxon is observed at the position x' in a given injection 
cycle (Fig. 4), the evolution of the amplitudes Cj of the 
measured system due to interaction with this fluxon is: 

1/2 

c j ■ ->^o{x' -x 3 )c /[Y>j\cjTp {x' -Xj)\ 2 } . (48) 

Qualitatively, and similarly to the conditional evolution 
in the transmission mode, the sequence of transforma- 
tions l|48|l describes "weak measurement". The system 
wavefunction is localized gradually in one of the states 
\j) with increasing number of the fluxon scattering events 
which lead to accumulation of information about In 
contrast to the transmission-detection mode, conditional 
evolution 148|) always remains pure under the simplifying 
assumptions adopted in this work for the JTL detector: 
coherent propagation of fluxons with fixed wavepacket 
ipo{x) and the same velocity. Fluctuations in the fluxon 
generator or finite dissipation in the JTL would create 
information losses also in this regime. 



V. NON-QUANTUM-LIMITED DETECTION 

Quantum-limited operation of the JTL detector dis- 
cussed in the preceding Sections requires, in essence, 
quantum-coherent dynamics of fluxons in the JTL. While 
this dynamics can be observed experimentally |2^, the 
task of realizing it is certainly very difficult. From this 
perspective, it is important that several attractive fea- 
tures of the JTL detector, e.g., large operating frequency 
and reduced parasitic dephasing during the time inter- 
vals between the fluxon scattering, remain even in the 
non-quantum-limited regime. Although the Josephson 
junctions of the JTL and those in the "external" parts of 
the JTL detector (generator and receiver) can give rise 
to dissipation and dephasing not related directly to mea- 
surement, in the JTL geometry (Fig. 1), parasitic dissi- 
pation is suppressed due to screening by the supercurrent 
flow in the JTL junctions [l4j . 

The dominant deviation from the quantum-limited de- 
tection should be associated then with the fluctuations 
in the fluxon motion. These fluctuations make the de- 
phasing factor v (I20|l due to scattering larger than the 
amount of information conveyed by the fluxon scattering. 
Some interesting measurement strategies are still possible 
with the JTL non-ideality of this type. The most natural 
example is the quantum non-demolition (QND) measure- 
ments of quantum coherent oscillations in a qubit mm 

which are designed to make the back-action dephasing by 
the detector irrelevant. Dynamics of the fluxon scatter- 
ing in the JTL detector makes it particularly suitable for 
the "kicked" version of the QND qubit measurements of 
a qubit ^3 or harmonic oscillator |4lj . 
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Consider, for instance, a qubit with the Hamiltonian 



H = -(A/2)a x , 



(49) 



which performs quantum coherent oscillations with fre- 
quency A. The qubit is coupled to the JTL through 
its a z operator, i.e. the states \j) used above, are the 
two eigenstates of a z . We assume that the Hamiltonian 
(|49|l already includes renormalization of parameters due 
to the qubit-detector coupling. If the qubit oscillations 
are weakly dephased at the rate 7 -C A (e.g., by resid- 
ual parasitic dissipation in the JTL detector), the time 
evolution of the qubit density matrix p during the time 
intervals between the successive fiuxon scattering events 
can be written in the ^-representation as follows: 



Pit) = \ 



1 



,—yt 



X 

iy 



-w 

-x 



r = —iAr , r — x + iy , 



(50) 



(51) 



where r(t = 0) = ±1 depending on whether the qubit 
starts at t = from the a z — 1 or er, — — 1 state. 



\<5=\> 



-1> 



Jt/A 



1/f 



t 

FIG. 5: Schematics of the QND fiuxon measurement of a 
qubit which suppresses the effect of back-action dephasing on 
the qubit oscillations. The fiuxon injection frequency / is 
matched to the qubit oscillation frequency A: / ~ A/ir, so 
that the individual acts of measurement are done when the 
qubit density matrix is nearly diagonal in the cr z basis, and 
the measurement back-action does not introduce dephasing 
in the oscillation dynamics. 

The fiuxon scattering at times t n = n/ f leads to partial 
suppression of the off-diagonal elements of p: 



y(t n + 0) = vy(t n - 0) . 



(52) 



However, if the fluxons are injected in the JTL with the 
time interval 1/f close to the half-period 7r/A of the 
qubit oscillations (Fig. 5), the qubit density matrix (|5*T|l is 
nearly diagonal at the moments of scattering, y(t n ) -C 1, 
and suppression (|52fl does not affect the qubit strongly. 
Such a QND measurement is possible with the JTL de- 
tector operating in any detection mode; to be specific we 
assume the transmission mode. The dependence of the 
fiuxon transmission probability on the qubit state can be 
written then as T + a z ST, where ST <C T in the "linear- 
response regime" . For quantum-limited detection, the 
linear-response condition ST <C T implies that v — > 1. 
In the non-quantum-limited case, the back-action can be 



stronger, and we take v to be arbitrary within the [0, 1] 
interval. 

If the frequency / is matched precisely to the qubit 
oscillations, / = A/-7T, the detector does not affect the 
qubit dynamics at all. If the mismatch is non-vanishing 
but small, 8 = A/ f — 71" <C 1, diagonal elements of p l|51|) 
evolve quasi-continuously even if suppression factor v is 
not close to 1. Equations 151|) and l|52|l give the following 
equation for this quasi-continuous evolution: 



x = -[(l + v)/(l-u)](fS 2 /2)x. 



(53) 



In the assumed linear-response regime, the qubit os- 
cillations manifest themselves as a peak in the spectral 
density Sy(w) of the voltage V across the JTL junc- 
tions. For / ~ A/-7T the oscillation peak in Sv{oj) is at 
zero frequency. Equations (f5T7|) and lt53|) describing the 
decay of correlations in the qubit dynamics in the a z ba- 
sis imply that the oscillation peak has Lorcntzian shape: 



S v (lo) = S + 



2Yp{ST)^l 



r 2 



and the oscillation linewidth T is: 

1+^(A-7T/) 2 



r = 7 - 



1 



2/ 



(54) 



(55) 



For the quantum- limited detection, v — > 1, and 
Eq. (|55fl reproduces previous results for the QND mea- 
surement | 40| , if one introduces the back-action dephas- 
ing rate Y ^ = f(l — v). In this case, Eq. (|55|) is valid 
for sufficiently small mismatch between the measurement 
and oscillation frequencies, \S\ <£L 1 — v. For larger 8, the 
oscillation peak in the detector output Sv(u>) moves to 
finite frequency A — nf and the QND nature of the mea- 
surement is lost 39] . Equation l|55|l shows also that in the 
limit of "projective" measurements v = 0, the broadening 
of the oscillation peak is weaker, and the peak remains 
at zero frequency for all reasonable values of the detun- 
ing parameter \8\ <C 1. Therefore, the stronger back- 
action of the JTL detector is advantageous for the QND 
measurements of coherent oscillations. The last remark 
is that although our discussion here assumed that the 
fiuxon arrival times are spaced exactly by 1/f, Eqs. I|54|l 
and H55|) should remain valid even in the presence of small 
fluctuations of the measurement times. These fluctuation 
can be described by taking into account that the detun- 
ing I S I cannot be made smaller that the relative linewidth 
of the fiuxon generator. 



VI. CONCLUSION 

We have suggested and analyzed a new ballistic detec- 
tor for the rapid read-out of flux qubits which can be 
implemented with the present-day SFQ fabrication tech- 
nology. The detector should combine quantum-limited 
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dynamics with time resolution that is better than the de- 
coherence time of typical flux qubits, making it possible 
to perform measurements on quantum-coherent qubits. 
To operate in the quantum limited regime, the detector 
requires junctions with sub- micron size, and its time reso- 
lution in this regime should reach the range of 0.3 ns. The 
detector is potentially useful even in the non-quantum- 
limited regime, when the time resolution can be improved 
further, and the detector can perform, for instance, the 
QND measurements of coherent oscillations. 
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